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In the high dimension (mean field) limit the susceptibility and the second moment correlation 
length of the Ising ferromagnet depend on temperature as x(^) — r ~ 1 an d £{T) = j 1 ^ 1 / 2 T - 1 / 2 
exactly over the entire temperature range above the critical temperature T c , with the scaling variable 
r = (T — T c )/T. For finite dimension ferromagnets temperature dependent effective exponents can 
be defined over all T using the same expressions. For the canonical two dimensional square lattice 
Ising ferromagnet it is shown that compact "extended scaling" expressions analogous to the high 
dimensional limit forms give accurate approximations to the true temperature dependencies, again 
over the entire temperature range from T c to infinity. Within this approach there is no cross-over 
temperature in finite dimensions above which mean-field-like behavior sets in. 

PACS numbers: 



I. INTRODUCTION 

The remarkable critical behavior at second order phase 
transitions has been intensively studied for many years. 
In the limit where the temperature T tends to the critical 
temperature T c , observables Q{T) diverge as 

Q{T)~Q\t\-i (1) 

where q is the critical exponent and t is the scaling vari- 
able; the textbook approaches to use t = (T — T c )/T c 
as the scaling variable, with small corrections, both ana- 
lytic and non-analytic, to the strict critical form as soon 
as T is not infinitesimally close to T c . It is widely consid- 
ered that there is only a narrow temperature range, the 
critical region, in the immediate vincinity of T c , where 
Eqs. ([!]) are valid and that they break down beyond this 
region. Thus the standard protocol for estimating criti- 
cal exponents from experimental or numerical data is to 
carry out analyses using Eqs. (|T|) with scaling variable 
i, together with finite size scaling (FSS) rules derived 
using these equations, over as narrow a range of temper- 
ature as possible around the critical point, introducing 
phenomcnological corrections to scaling if the data are 
accurate enough. 

In the discussion that follows we will consider the sec- 
ond moment correlation length £, sm (T) (not the "true" 
correlation lengthP £ trlle ) which we will refer to as £(T), 
and the reduced susceptibility x(T) (see ^ for the defi- 
nitions). It should be kept in mind that the standard 
critical observable, the reduced susceptibility x(T), is 
by definition related to the thermodynamic susceptibility 
Xth = (dM /(IH)h^o through the temperature dependent 
normalization x(T) = T\th(T). 

From arguments based on high temperature series ex- 
pansions (HTSE)0 an " extended scaling" formulation was 
introduced^ in which to leading order the reduced sus- 
ceptibility is written x(^) ~ [ T ] 7 an d the second mo- 
ment correlation length is written £(T) ~ r~ 1 / 2 [r]~", 
with the scaling variable r = (T — T c ) /T^. Here we will 



discuss two extreme canonical systems : the high dimen- 
sional near neighbor Ising ferromagnet on a hypercubic 
lattice, and the two dimensional Ising ferromagnet on a 
square lattice. In the high dimensional (mean field) case 
the critical behavior parameterized this way is strictly 
exact with temperature independent effective critical ex- 
ponents over the entire temperature range from T c to 
infinity. In the square lattice case the approach leads to 
compact accurate approximations over the whole range. 
In the temperature range close to T c the expression for 
£(T) can be linked to Renormalization Group Theory 
(RGT) analytic corrections. Using this protocol, exper- 
imental or simulation data can be usefully analysed in 
terms of the critical behavior over a wide temperature 
range and not only in the very close neighborhood of T c . 

II. HIGH DIMENSION ISING MODEL 

It is instructive to first consider the HTSE for the high 
dimension Ising ferromagnet in the infinite dimension (or 
mean field) limit. For near neighbor interaction spin 
1 /2 Ising ferromagnets on [hyper] cubic lattices the HTSE 
for the reduced susceptibility-^ x(P) an d for the second 
moment of the correlation function /J-2(P) are given ex- 
actly by infinite series of terms in powers of tanh{(3) with 
(3 = l/T: 

x (/3) = 1 + (z)tanh{/3) + (z 2 - z)tanh(P) 2 

+(z 3 - 2z 2 + z)tanh((3) 3 + ■■■ (2) 

and 

M2 (/3) = ( z )tanh{(3) + (2z 2 - z)tanh(f3) 2 

+ (3z 3 - 2z 2 + z)tanh{(3f + ■■■ (3) 

where z is the number of neighbors (z = 2d for hyper- 
cubes in dimension d). 

The second moment correlation length is defined by 

m = [H2{P)/ZX{P)] 1/2 (4) 
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When z — » oo the (ztanh(P)) n contribution will domi- 
nate at each n, and (3 C ~ l/z will decrease so tanh(fl) — > 
/? in all the paramagnetic regime. Hence 

x(/?) = l + z/3+(z/3) 2 + (z/3) 3 + --- 

= (l-z/3)- 1 (5) 

and 

M2 (/3) = zx(/?)£(/?) 2 
= + 2(zp) + 3(z(3) 2 + A(z(3) 3 + • • •] 

= 2/3(1 - z/?)- 2 (6) 

i.e. 

x(/3) = (i-/3//3 c )^ 1 -[(r-r c )/r]- 1 (7) 

and 

£(/3) = ^{l-P/Pc)- 1 ' 2 = T-WftT-Tc)/!]- 1 ' 3 (8) 

exactly for all /3 less than (3 C = 1/z, i.e. all T greater 
than T c = z. 

Using as the critical variable r = (1 — (3/(3 c ) instead 
of i and introducing the prefactor /3 1 / 2 in the expression 
for £(/3), the critical regime as defined by the scaling ex- 
pressions : 

x(/3) = r-T (9) 

and 

m/p 1 ' 2 = r-" (io) 

extends rigorously from T c to infinite T, with /3 C = l/z 
and temperature independent mean field exponents 7=1 
and v = 1/2. 

This statement can be reformulated : in the "ideal" 
mean field limit one can expect correction terms to be in- 
existant if the critical variable and observables are chosen 
correctly. The two observables which we have discussed, 
the reduced susceptibility x(/3) (i- e - Xth(P)/{3) and the 
"reduced" second moment correlation length £(/?) / /3 1 / 2 , 
show exact critical power law behaviors (as r -7 and as 
t~ v respectively) for all T > T c . We can surmise that 
in finite dimensions there will be correction terms but 
that the same critical variable and normalized observ- 
ables will remain the most appropriate for expressing the 
critical behavior over a wide temperature range. 



III. SQUARE LATTICE ISING MODEL 

In finite dimensions the extreme simplicity of the mean 
field case will be lost, but because of the generic struc- 
ture of the HTSE shown above the general form of equa- 
tions Eqs. (pi and ( |10[ ) (including the non-critical normal- 
ization /3 1 "in C(/3j//3 1/2 ) can be expected to be rather 



robust 4 . The exact finite dimension leading critical be- 
haviors for [hyper] cubic ferromagnets when (3 — » f3 c can 
be written 

X (f3) - C x (t)-< (11) 

and 

m -> [Q/4 1/2 ]/3 1/2 (r)- v (12) 

1 /2 

with critical amplitudes C x , C^//3 C ■ For spin 1/2 the 
infinite temperature (3 — > values are x(/3) ~^ 1 and 
£(/?) /3 1/2 for all d. Extended scaling expressions can 
be written which are analogous to the infinite dimension 
form, linking the critical limit with the trivial high tem- 
perature fixed point limit!^. As a first step a minimal 

modification must be made in order to allow for the fact 

l li 

that in finite dimensions C x and C^/(3 C are not equal 
to exactly 1. We write : 

X*W = C x [t]-<[1 + t(1-C x )/C x ] (13) 

and 

f(/3) =/3 1/2 [(Q//3 c 1 / 2 )(r)-][l + r(/3 c 1 / 2 -Q)/Q] (14) 

These expressions, which depend only on the critical 
parameters /3 C >7) v> C x and Cg, are exact by construction 
in both the critical and the high temperature limits, and 
they provide compact approximate expressions for the 
behavior over the whole range in between. It has been 
demonstrated that for the standard thre e dimensional 
Ising, XY and Heisenberg ferromagnets^' extended scal- 
ing expressions (defined entirely through the critical pa- 
rameters appropriate for each particular case) agree with 
the true x(T) and £(T) to within better than about 1% 
over the entire range of temperature from T c to infinity. 

Here we will consider in more detail the particular case 
of the canonical two dimension square lattice Ising fer- 
romagnet, for which the critical temperature, the crit- 
ical exponents, and a number of other properties are 
known exactly from the original work by Onsager and 
others^ and from more recent conformal field theory^. 
In the square lattice there are no irrelevant operators^ 
except those due to the lattice breaking of rotational 
symmctr ySa which can be ignored for present purposes. 
There are no exact analytic expressions for the suscep- 
tibility or the second moment correlation length but 
many ter ms of th e high temperature series expansions 
are knowr l 11 ' 12 ' 1 ^ so x(/?) and £(/3) can be calculated 
to high precision over the entire range of (3, from (3 C to 
zero (i.e. from T c to infinity). We have taken advan- 
tage of the exact knowledge of the critical temperature 
tanh([3 c ) = s/2 - 1 (i.e. (3 C = 0.440686794 • • •) and of 
the critical exponents 7 = 7/4 and v = 1 to form biased 
Pade approximants of x and £ 2 , using up to 48 HTSE 
coefficients. Using first-order inhomogencous differential 
approximants we would obtain an equivalent accuracy. 
Thus C x = 0.962581 • • ■, and C ( = 0.854221175. (The 
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ratio of the exact "true" correlation length critical am- 
plitude to the second moment correlation length critical 
amplitude is 1.000402-31). The extended scaling estimates 
X*(/3) and £*(/?) (Eqs. (p) and (fl4)) can be written 



down directly using these values. 

Fisher and BurforcP some 40 years ago introduced a 
non-critical factor in the expression for £(/3) for Ising fer- 
romagnets ; they already noted that in the mean field 
limit the prefactor would be equal to /3 1 / 2 as confirmed 
by the high dimension limit discussion above. For Ising 
ferromagnets in dimensions 2 and 3 they introduced an 
"effective range of direct interaction" parameter ri(T), 
which they defined through 



MT)/£(T)) 
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1/X(T) 



(15) 



The temperature variation of r\ (T) is mainly due to the 
non-critical prefactor in the expression for the second mo- 
ment correlation length Fisher and Burford did not 
give an explicit expression for r\(T) in finite dimensions 
but they calculated it numerically over a wide range of 
temperature from the HTSE terms known at the time 
for five different Ising ferromagnets in dimensions 2 and 
3 (their Figure 6) . We can recalculate the Fishcr-Burford 
square lattice >"i(T) as defined above from the high pre- 
cision square lattice Pade approximants values and also 
from the extended scaling expressions. In Figure 1 we 
compare over a very wide temperature range 



^i(r)=C P (r)/(x P (r)) 1/(2 -" ) 



(16) 



obtained directly from the high precision values 
£, p (t) , x p (t) with rl(r) estimated from the £*(t),x*(t) 
of the extended scaling expressions. On the scale of this 
plot the two sets of points are almost indistinguishable, 
showing that the extended scaling expressions are accu- 
rate approximations to the exact behavior of the observ- 
ables. 

As a more critical test of the extended scaling, in 
figures 2 and 3 we present the ratios £*( r )/£p( T ) an d 
X* ( T ) / Xp( T ) where again £*(t), X*(t) are the extended 



scaling estimates from equations Eq. ( 14 1 and Eq. ( 13 1 



and £p( r ); Xp( T ) arc the high precision calculated values. 
It can be seen already that without further correction 
factors the extended scaling values represent the high 
precision data for x(r) and for £(r) to better than 1% 
over the entire temperature range from r = to r = 1. 
The figures also include £*(t)/£ p (t) and X*( T )/Xp( T ) ra- 
tios where the expressions for £*(r) and x*( r ) include 
second correction terms from Eqs. [22] and [23] as will be 
discussed later. 

An alternative manner in which to present the data is 
to express temperature dependencies of x( T ) an d £,(t) 
in terms of tem perature dependent "effective" critical 
exponent d 3 * 6 * 16 l In the spirit of the previous discussion, 
quite generally for any ferromagnet the temperature de- 
pendent effective exponents can be rigorously defined as 
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FIG. 1: The Fisher-Burford "effective range of direct inter- 
action" parameter r\(T) calculated from the high precision 
square lattice Pade approximants values (black squares) and 
from the extended scaling expressions (red circles), as func- 
tions of t = 1 — /3/f3 c . 
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FIG. 2: The ratio £* (t)/£ p (t) where the £*(t) are the ex- 
tended scaling estimates of the correlation length (red circles 
one correction term, blue triangles two correction terms) and 
£ p (t) the high precision data, as a function of r = 1 — @//3 c . 



u eff (r) = dlog(t(P)/P 1/2 )/dlog(T) (18) 



and 



leff{r) = dlog{x{r)) / dlog{r) 



(17) 



Tfc//G9) = 2 - dlog( x ((3))/dlog(m/P 1/2 ) (19) 

The definition of J e ff(T~) is the same as in references^!, 
but because of the prefactor /3 1 / 2 normalizing £(/?) the 
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FIG. 3: The ratio X* ( T ) I Xv( T ) where the X*( T ) are the ex- 
tended scaling estimates of the reduced susceptibility (red 
circles one correction term, blue triangles two correction 
terms) and Xp( T ) the high precision data, as a function of 

T = 1 - f3/f3 c . 



FIG. 4: The effective exponent Jeff calculated from the high 
precision Xp( T ) data (black squares) and from the extended 
scaling estimates ( red circles one correction term, blue trian- 
gles two correction terms), as functions of r = 1 — /3//3 c - 



other two definitions are not standard. This prefactor 
is essential to insure sensible high temperature limits 
in Eq. [15] and Eq. [19] We can note that for the first 
two parameters an explicit choice must be made for /3 C , 
while rj e ff{(3) (obviously linked to the two others) can 
be calculated from x(P) and data sets without any 



a priori knowledge of or estimate for j3 c as Eq. [19] does 
not involve r. The effective exponents tend to the crit- 
ical exponent values at f3 c , and from HTSE there are 
simple exact results for the high temperature limits : 
7e//(0) = 2d/3 c , v e ff(Q) = d(3 c and r? e //(0) = in all 
dimensions for hypercubic lattices. In the square lattice 
case 7e//(0) = 1.7627- • • and v eff (Q) = 0.8814- • -. 

Figures 4, 5 and 6 show J e ff(r), v e f f (r), and r} e ff(r), 
comparing the high precision values from HTSE with the 
extended scaling values, again with one or two correction 
terms. 

A number of remarks can be made. First, looking only 
at the high precision data it can be seen that the effec- 
tive exponents defined through Eqs. (17 1, (18 1 and (19) 



change smoothly and gradually with temperature over 
the whole range of temperature from the critical temper- 
ature to infinity. For the square lattice 7e//0") changes 
little with temperature while v e jf(r) varies rather more. 
In the 2d case ry e //(r) must necessarily change quite 
strongly, from the critical value 0.25 to the infinite tem- 
perature value of 0. In dimension 3 the absolute value of 
the change is much weaker. 

Secondly, presenting the data in this way is a very sen- 
sitive test of the extended scaling expressions; it can be 
seen that with one correction term agreement with the 



high precision data is to within better than about 5% for 
each of the effective exponents over the entire temper- 
ature range. With two correction terms the agreement 
is considerably improved, particularly for v e ff- From 
the general argument given above, in other ferromagnets 
one should expect qualitatively similar behavior, with 
marginal complications due to weak irrelevant operators 
terms. 

A widely held concept is that one should expect a 
"cross-over" to mean- field- like effective exponents (7 = 
1,1/ = 1/2,1] — 0) outside some putative critical region 
in the neighborhood of (3 C . With the present definitions 
of the temperature dependent effective exponents there 
is clearly no such crossover; J e ff and v e ff do not tend to 
mean-field- like values, although rj e ff does head gradually 
to the mean field value of 0. 

(Even denning effective exponents in the traditional 
way using t as the scaling variable and no prefactor would 
lead to limiting high temperature values Jeff (tec) = 0, 
Veff (too) = 1/2, and 7? e //(*oo) = 00. Jeff (too) and 
Veff (too) are completely different from the mean field ex- 
ponents) . 



IV. ANALYTIC CORRECTIONS 

The general formalism for corrections to scaling 
(both analytic and non-analytic) within the Renormal- 
ization Group Theory (RGT) is well established (see 
references! 1211 ^. In the square lattice Ising model for 
present purposes it can be c onsidered that there are no 
non-analytic correction s' 9 * 10 !. This means that one can 



5 



1.00-1 




1 1 1 1 1 1 1 1 1 (— 

0.0 0.2 0.4 0.6 0.8 1.0 

X 



FIG. 5: The effective exponent v & /f calculated from the high 
precision £ p (t) data (black squares) and from the extended 
scaling estimates ( red circles one correction term, blue trian- 
gles two correction terms), as functions of r = 1 — /3//3 c . 
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FIG. 6: The effective exponent rj e ff calculated from the high 
precision Xp( T ) an< i £p( T ) data (black squares) and from the 
extended scaling estimates (red circles one correction term, 
blue triangles two correction terms), as functions of r = 1 — 

P/0C- 

write for instance the ratio of the inverse of the exact 
"true" correlation length l/£t™ e (r) = ln(coth(T)) — 2r 
to the ideal pure critical power law value r/C| as an an- 
alytic correction factor g(r) which is a smooth function 
having an asymptotic Taylor expansion form, i.e. 

C|/6™ e (r) = Tg(r) = r[l+aiT+a 2 r 2 +a 3 r 3 + - • •] (20) 



The angular dependence of the exact leading coefficients 
of this "large distance 2-point correlation length" ana- 
lytic correction series has been calculateo™^. If g(r) 
is truncated at low order it only gives a useful represen- 
tation of £,true(i~) very close to T c . 

For £ (it is important to keep in mind that this second 
moment correlation length is not the "true" correlation 
length) we choose to define a correction factor g*(r) in 
the extended scaling form : 

W = ((l-r)- 1/2 rYW (21) 

as (P/fic) 1 / 2 = (1 — t) 1 / 2 . Empirical fits can be made to 
the high precision £ p (t) square lattice data. We assume 
a Taylor series for <7*(t), keeping only two terms (up to 
order r 2 ) with the restriction that <?*(1) = Cj. We find 
that 

5 *(t) = 1 - 0.19055t + 0.04476-r 2 (22) 

provides a fit accurate to within less than 0.02% over the 
entire temperature range from t — to 1, thus includ- 
ing both the "critical region" and the high temperature 
region. If the standard g{r) correction factor had been 
used the coefficient for the term in r would have become 
0.5 - 0.19055 = 0.30945. This value is very close to the 
value of the equivalent coefficient in the g{r) series for 
the "true" correlation length, which is 0.3116. 

The large coefficient of the r term in g*(r) can 
be understood as correcting for the difference between 
C* ? //3 C 1/2 = 0.854221 and 1. The small second coefficient 
corresponds to a further adjustment which leads to a very 
significant improvement in the fit. It is remarkable that 
excellent precision over the entire temperature range is 
obtained with only two correction terms, validating the 
use of the g*(r) parametrization rather than the g{r) 
form for ^ sm (r). We can note that when g*(r) is used 
rather than g(r) there are no correction terms to £ sm (T) 
in the mean field limit. 

The square lattice correction terms for x(r) have been 
carefully studied-^. The coefficients for six leading cor- 
rection terms are known of which three are exact with 
additional terms over and above Taylor series terms. An 
empirical correction factor for x( T ) with two correction 
terms gives 

x (t) = C x r- J/4 (l + 0.0779r - 0.03903t 2 ) (23) 

The leading term in Eq. [23] has a coefficient essentially 
identical to the exact value (0.0779032 ■ • •) and the phe- 
nomenological inclusion of a single further term provides 
an overall fit up to infinite temperature which is accurate 
to better than 0.1%. 



V. FINITE SIZE SCALING 

An immediate practical consequence that follows from 
the discussion above concerns the extraction of critical 
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parameters from numerical studies of systems other than 
these canonical ferromagnets. If numerical data have 
been obtained over a wide temperature range and not 
only in the region very close to T c , direct plots of the 
effective exponents as defined through Eqs. (17 1 to (19) 



for the largest sample sizes available can be extrapolated 
to estimate critical exponents. The temperature depen- 
dence of the curves, including the temperature range well 
above T c , can then give useful indications concerning crit- 
ical exponents and corrections. Finite size scaling analy- 
ses should be made using appropriate expressions derived 
from the leading extended scaling form. The widely used 
finite size scaling relation 



Q(L, T) ~ F[L}I V (T - T ( 



(24) 



is derived from the Fisher finite size scaling ansatz 
Q(L,T) = F[L/£(T)] on the assumption that the cor- 
relation length behaves as £(T) ~ t~ u ; as we have seen 
the latter is always a very poor approximation for £(T) 
except extremely close to T c . Using the Fisher ansatz 
together with the extended scaling rule ~ (3 1 I 2 t~ v 

leads to finite size scaling expressions 4 ^ which should 
remain much better approximations over a considerably 
wider temperature range. If T c is known to reasonable 
precision such finite size scaling analyses from numerical 
data taken over a wide temperature range should give re- 
liable and unbiased estimates for the critical exponents. 
Ideally an allowance for the correction factor g* (r) of the 
preceding section should also be included in the analysis 
but this would require very high quality data. 

It is important that for each particular system the ap- 
propriate extended scaling form should be used. For in- 
stance, in spin glasses with symmetrical interaction dis- 
tributions the right scaling variable i d 4 * 17 ! 

tsg = i - om) 2 . 



VI. CONCLUSION 

One should expect "ideal" critical behavior in the mean 
field limit ferromagnct, meaning that if the scaling vari- 
able and the normalizations of the observables are chosen 
appropriately all observables should show pure critical 
power law behavior over the entire temperature range 
above T c . In the high dimensional limit both the re- 
duced susceptibility x( T ) an d the "reduced" second mo- 
ment correlation length £ (t) //3 1 / 2 indeed show exact crit- 
ical power law behaviors (as r f and as t~ v respectively 
with temperature independent mean field exponents) for 
all T > T c , validating the use of the scaling variable 
r = (1 — fl/Pc) and the normalization for £(r) through 



the (3 1 / 2 non-critical prefactor. In other words, in terms 
of this scaling variable and these observables behavior is 
" always critical" for the mean field system over the whole 
temperature range. 

For systems in finite dimensions one can no longer ex- 
pect ideal behavior to be followed exactly at all tem- 
peratures, and there will always be deviations from the 
pure critical power laws as soon as r is finite; these de- 
viations can be expressed in terms of correction factors. 
However in order to reduce the importance of the neces- 
sary corrections (and to have sensible high temperature 
limits) it is judicious to base the choice of scaling vari- 
able and scaling expressions on those appropriate for the 
mean field limit. In agreement with arguments from the 
general form of the HT series expansions, this leads to ex- 
tended scaling expressions with a single correction term 
which are X *(P) = C x [l - + r(l - C x )/C x ] 

for the reduced susceptibility and = /5 1 / 2 [Q(1 — 

/3//3 c ) -i/ ][1+t(1— C^/Cy for the second moment correla- 
tion length. The adjustment terms have been introduced 
in order that the expressions tend to their known exact 
high temperature limits. For the canonical square lattice 
Ising ferromagnet with its known critical temperature, 
exponents, and critical amplitudes C' x and C^, the ex- 
pressions x*{fl) an d £*(/?) are compact approximations 
to the exact behavior which are accurate to within about 
1% over the entire temperature range. If the correction 
series is truncated after only one more term (so that it 
takes the form [1 + a\T + CL2T 2 )) the precision improves 
to better than 0.1% and better than 0.02% for x(/3) and 
£(/?) respectively. 

The temperature dependence of the observables can be 
expressed in terms of strictly defined temperature depen- 
dent effective exponents, Jeff( T )> v zff(j) an d v eff(j), 
which vary smoothly and weakly with temperature. J e ff 
and v e ff do not tend to mean field values at high tem- 
perature. In the square lattice system the two term cor- 
rection factor for £*(/?) can be compared with the critical 
analytic correction factor for the "large distance 2-point" 
correlation length, for which the exact leading terms in 
the Taylor expansion have been discussed in detail in the 
RGT formalism^. 

Within the extended scaling approach there is no cross- 
over temperature for finite dimension systems above 
which mcan-ficld-likc behavior sets in. 

The present conclusions confirm those already given 
following analyses of data from other ferromagnets and 
spin glasses^. For practical purposes these conclusions 
should be taken into account when analysing experimen- 
tal results or when extracting critical parameters from 
finite size scaling analyses on numerical simulation data. 
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